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ABSTRACT  Detailed  discussions  are  made  of  principles  and 
methods  associated  with  orbital  determinations  for  low  orbit 
spacecraft  making  use  of  GPS  systems.  Analysis  is  done  of 
problems  which  should  be  considered  with  regard  to  low  orbit 
spacecraft  during  selection  of  satellites.  In  conjunction  with 
this,  satellite  selection  tests  are  carried  out  on  the  basis  of 
navigation  satellite  as  well  as  user  models.  Discussions  are 
made  of  making  use  of  least  square  methods  to  make  calculations 
to  determine  orbits.  Detailed  analysis  is  done  of  Kalman  filter 
algorithms.  With  a  view  to  low  orbit  spacecraft,  a  nine  state 
Kalman  wave  filter  model  is  put  forward.  In  conjunction  with 
this,  simulation  analog  calculations  are  carried  out  on  this 
foundation.  Besides  this,  on  the  basis  of  special  user 
requirements  with  regard  to  electric  power  source  energies  and 
reliability,  one  type  of  relatively  practical  orbit  determination 
method  is  put  forward.  Simulation  results  clearly  show  that  the 
method  put  forward  satisfies  user  requirements  in  terms  of 
accuracy  and  is  appropriate  for  actual  engineering  applications. 

SUBJECT  TERMS  Global  positioning  system  Least  square 
method  Kalman  filtering  Satellite  algorithm 


1 


The  global  positioning  system  (GPS)  is  a  second  generation 
satellite  navigation  system  developed  on  the  foundation  of  the 
TRANSIT  satellite  navigation  system.  At  the  present  time, 
installation  work  on  the  whole  24  satellite  navigation 
constellation  is  already  completed.  It  is  possible  to  realize 
contiuous,  real  time,  all  weather  three  dimensional  positioning 
on  a  world  scale.  There  are  broad  applications  in  such  realms  as 
land  communications,  geodetic  surveys,  map  making,  maritime 
navigation,  aviation,  and  so  on.  The  prospects  for  applications 
in  the  field  of  space  are  also  very  wide. 

At  the  present  time,  there  are  basically  the  three  types  of 
methods  below  for  determining  spacecraft  orbits  [1]:  ground 
equipment  orbit  determination  method,  tracking  and  data  relay 
satellite  system  (TDRSS)  orbit  determination  method,  and  the  GPS 
orbit  determination  method.  Looking  at  the  situation  of  China, 
due  to  limitations  associated  with  geographical  breadth  and 
economic  conditions,  making  use  of  the  first  two  types  of  methods 
has  definite  limitations  in  both  cases.  As  a  result,  making  use 
of  GPS  autonomous  orbit  determinations  is  one  type  of  very  good 
choice.  This  article  discusses  methods  associated  with  making 
use  of  Kalman  filter  algorithms  and  least  square  methods  to 
calculate  the  orbits  of  low  orbit  spacecraft.  In  conjunction 
with  this,  imitative  simulations  are  carried  out. 

1  ALGORITHMS  FOR  SELECTING  SATELLITES 

As  far  as  making  use  of  GPS  to  carry  out  navigational 
positioning  is  concerned,  users  normally  want  to  select  four  or 
more  navigational  satellites.  The  relative  geometrical  positions 
of  these  four  navigational  satellites  has  very  great  influences 
on  navigational  positioning  accuracy.  When  ranging  errors  are 
fair  si^ed,  positioning  errors  produced  by  different  geometrical 
relationships  between  users  and  these  four  satellites  are  also 
different.  As  a  result,  selection  should  be  made  of  satellite 
constellations  associated  with  optimum  geometrical  positions  in 
order  to  obtain  the  best  positioning  accuracy. 

Up  to  the  present  time,  GPS  system  satellite  constellation 
selection  methods  which  relevent  references  have  introduced  can 
be  divided  into  the  four  types  below  [2]:  optimum  geometrical/23 
precision  factor  method  [3],  maximum  vector  terminal  tetrahedron 
volume  method  [4],  maximum  orthogonal  projection  method  [5], 
combination  method  [6].  These  four  types  of  methods  each  have 
their  respective  characteristics.  With  regard  to  ground  or  low 
altitude  targets  (aircraft,  and  so  on) ,  all  of  them  are  capable 
of  obtaining  relatively  good  results.  However,  speaking  in 
regard  to  low  orbit  spacecraft  several  hundred  thousand  meters 
off  the  ground — due  to  limitations  in  the  radiation  angles  of 
antennas  on  GPS  navigation  satellites  (<  ±  14.3°) — user  receiver 
antenna  directional  patterns  are  not  complete  either  (generally 
smaller  than  ±80°).  Besides  this,  there  are  also  such  problems 
as  the  blocking  of  solar  energy  battery  sail  panels,  and  so  on. 
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As  a  result,  the  selection  of  navigation  satellites  still  has 
some  special  spots  to  it.  Below,  a  usability  analysis  is  carried 
out. 


Fig.l  "Usable"  Navigation  Satellite  Schematic 

GPS  navigation  satellites  are  20182km  away  from  the  ground. 
Due  to  limitations  on  their  antenna  radiation  angles,  and  user 
receiver  antenna  directional  pattern  angles  smaller  than  ±80°,  as 
a  result,  when  low  orbit  spacecraft  are  flying  at  a  certain 
altitude  in  space,  there  is  a  possibility  to  create  "visible"  and 
not  "usable"  situations.  That  is,  on  the  basis  of  receiver 
antenna  directional  patterns,  navigation  satellites  are  indeed 
located  within  the  receiving  range,  but  the  transmited  beams  of 
navigation  satellites  are  still  not  able  to  cover  user  receiving 
antennas.  Because  of  this,  the  navigation  satellites  in  question 
are  then  not  "usable".  Navigation  satellites  which  can  be  used 
to  supply  user  positioning  should  possess  two  conditions — they 
must  be  "visible" ,  and  they  must  also  be  "usable" .  As  is  shown 
in  Fig.l,  when  p  <  80°,  navigation  satellites  are  visible.  When 
a  <  143°,  navigation  satellites  are  usable.  When  the  geometrical 
positions  of  navigation  satellites  and  users  simultaneously 
satisfy  the  two  conditions  above,  can  just  be  used  in 
positioning.  In  Fig.l,  satellite  S  is  usable  and  visible  with 
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respect  to  u.  However,  it  is  visible  but  not  usable  with  regard 
to  u'. 

After  precisely  determining  usable  satellites,  it  is  then 
possible  to  opt  for  the  use  of  various  types  of  methods  to 
calculate  GDOP  values  in  order  to  precisely  specify  four 
navigation  satellites  used  for  navigational  positioning. 

2  DATA  PREPROCESSING 

GPS  systems  are  passive  ranging  systems.  Users  go  through 
measurements  of  their  own  relative  to  already  known  positions, 
speeds,  and  clock  differences  associated  with  pseudo  ranges  and 
pseudo  range  rates  of  the  four  navigation  satellites  to  precisely 
determine  three  dimensional  coordinates,  speeds,  and  clock 
differences  of  user  clocks  relative  to  GPS  time.  On  the  basis  of 
analysis,  one  has: 


In  this,  ap  is  positioning  error; 

a  is  the  mean  square  deviation  associated  with  pseudo 
range  measurements; 

GDOP  is  the  geometrical  error  factor. 

It  is  possible  to  see  that  positioning  errors  are  dependent 
on  geometrical  error  factors  as  well  as  the  precision  of  pseudo 
range  measurements.  In  order  to  increase  the  precision  of 
ranging,  on  the  one  hand — in  equipment — it  is  possible  to  opt  for 
the  use  of  various  types  of  methods.  On  the  other  hand — in  terms 
of  data  processing — it  is  possible,  by  contrast,  to  opt  for  the 
use  of  data  smoothing  methods  in  order  to  reduce  ranging  errors 
carried  along  with  random  noise. 

With  respect  to  internal  reciever  noise  as  well  as  exterior 
interference,  pseudo  ranges  obtained  in  measurements  can  appear 
strange,  that  is,  "wild  values".  If  preprocessing  is  not  gone 
through  and  direct  use  is  made  in  positioning  calculations, 
relatively  large  errors  will  be  produced.  This  article  opts  for 
the  use  of  conversion  methods  to  carry  out  preprocessing  with 
regard  to  measurement  values. 

First  of  all,  a  conversion  time  period  interval  (normally 
Is)  is  selected.  Each  pseudo  range  measurement  value  Sn(tj) 
within  this  time  interval  is  taken  (in  this,  tH  <  tj  <  tk  ,  tk  is 
the  interval  end  point)  and  converted  to  interval  end  /24 


In  this  k  —  carrier  wave  wave  length; 

Aijrn  (tj,tk)  —  time  interval  carrier  wave  phase  increment 

from  tj  to  tk  . 
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After  that,  the  mean  pseudo  range  value  is  calculated  out 
for  instant  tk 

=  (3) 


As  far  as  application  of  this  type  of  method  is  concerned,  it  is 
possible — to  a  very  great  extent — to  restrain  the  influences  of 
abnormal  values,  thereby  increasing  positioning  calculation 
accuracy. 


3  USE  OF  LEAST  SQUARE  METHODS  IN  ORBIT  DETERMINATION 


Processing  of  observation  data  is  one  important  link 
associated  with  measurement  operations.  Traditional  orbit 
determination  methods  make  use  of  least  square  principles.  After 
waiting  for  completion  of  observation  sampling,  use  is  made  of 
batch  processing  methods  to  calculate  the  orbit  associated  with  a 
certain  epoch  instant.  This  is  also  nothing  else  than  the 
orbital  improvement  which  is  often  spoken  of. 


psuedo  range 


r’  =r,  +  c  •  b 

=  [(x  —  x.)2  +  (y  —  y.y  +  (z  —  x,)2]1'2  +  c  •  b 

■ .  (x  —  x,)  (x  —  x,)  +  (y  —  y,)(.y  —  y,)  +  (z  —  z,)(z  —  z,) 
r ■  [(x  -  x,)2  +  (y  -  y,y  +  (z  -  z.yj'1 

pseudo  range 
rate 


(4) 


In  this,  x,y,z  are  three  dimensional  user  coordinates; 

x,y,z  are  three  dimensional  user  velocities; 
xs,ys,zs  are  three  dimensional  positions  of  navigation 

satellites; 

x.*y.'Z.  are  three  dimensional  velocities  of  navigation 
satellites; 

b  is  user  clock  deviation; 
f  is  user  frequency  deviation; 
c  is  the  speed  of  light. 

In  a  certain  estimated  user  state,  a  Taylor  expansion  is  carried 
out 
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r-  =/  (x,  y,  2,  x,  y ,  z,  b,  f ) 


=  r„+r*  •  Ai+o  (Ax2) 


*  Ay+^  •  Ax+f 

•  A£+|f  •  A/+o  (Ax2) 


A^+^T  *  A:v+ig-  ‘  Az+ 


In  these  equations,  o(Ax2)  is  a  high  order  term. 

After  omiting  high  order  terms,  the  matrix. is  written  in  the 
form  of  /25 


D{ 

Di 


Lr;  -  D.  J 


3x  dy 


£L  *L 

3x  dy 


L, 


3r i 

dz 


a rj_ 
dz 


dx 


dx 

A, 


?±±L 

dy  dz  3b  df 


a£  drf_  drf_  dr^ 
dy  dz  db  df 


'Ax' 
Ay 
Az 
Ax 
A y 
Az 
A  b 
AJ¬ 
AX 


(5) 


In  this,  n  is  an  integer  associated  with  pseudo  range; 

rio  =  Di  is  the  range  value  from  the  user  to  the  ith 

navigation  satellite  (i  =  l,2,...,n)  for  instant 
k; 

Di  =  [ (xsi  -  x)2  +  (ysi  -  y)  +  (zsi  -  z)2]l/2  ; 

x,  y,  z  are  user  coordinates  calculated  for 
instant  k; 
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In  this,  n  is  an  integer  associated  with  pseudo  speed; 

ru  =  A  which  is  the  relative  motion  velocity 

calculated  for  instant  k  between  the  user  and 
the  ith  navigation  satellite; 
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_  (x„  —  x)( x„  —  x)  +  (>v  —  ;y)(;y„  —  y)  4-  (s„  —  z)(z„  —  g) 

[(xw  -  x)2  +  (y,  -  >02  +  <z„  -  z)2]"2  ? 

x,y,z  are  three  dimensional  user  coordinates;  /26 
x,^,z  are  three  dimensional  user  velocities; 

=  - - —  +  [(x  —  x„)(x  —  x*)  +  (y  —  y „)(y  —  y„)  +  (z  —  z.)(z  —  z„) 

dx  r,  or 

^  4.  [(x  —  x„)(x  —  x.)  +  Cy  —  jy.-Xy  —  y.)  +  (z  —  z„)(z  —  z„)  v1; 

3y  r,  .  dy 

=  - - —  +  [(x  —  x„)(x  —  x„)  +  (>  —  —  S’.)  +  (z  —  z.)(z  —  z.) 

dz  r,  <» 

drj  _  x  —  x„  an*  _  >  —  31..  a^r  _  z  —  z„  a-;  ajv  =  , 

r,  ’  dy  r,  3z  r,  ’  c&  ’3/ 

On  the  basis  of  formula  (4)  and  formula  (6) ,  one  also  has 

an  _  x„-  -—  x  an_  - 
3x  r, 

ar,  _  >»„  ~  3>  _  3rJ_ 
dy  r,  dy 

dr l  _  z„  —  z  _ 

at  n  2s 


Merging  the  two  matrices  into  one,  one  obtains 


•  AX  =  A  •  AX 


(8) 


In  the  equations,  L  is  the  observation  matrix;  A  is  the  gradient 
matrix;  A  X  is  the  state  vector  correction  amount. 
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From  matrix  theory,  it  is  possible  to  know  that  the  equation 
A  ,  A  X  =  L’s  least  square  solution  is 

A  X  =  (ATA)-IAT  .  L 

From  this,  it  is  possible — on  the  basis  of  the  values  of 
each  iteration  of  measurements — to  calculate  user  state  deviation 
corrections  to  carry  out  calibrations  with  regard  to  the  next 
state,  Xk+1  =  Xk  +  AX.  After  that,  new  calculations  are  done  of  Di 
and  A  ,  making  use  of  new  observation  quantities  to 
calculate  L  matrices,  repeating  the  corrections  above.  In  this 
way,  after  going  through  a  certain  number  of  iterative 
substitutions,  it  is  then  possible  to  precisely  determine  user 
positions,,  speeds,  and  other  condition  values.  In  the 
calculations  discussed  above,  after  each  iterative  substitution, 
new  calculations  should  be  done  of  coefficient  matices  A  and 
observation  matrices  L. 


4  KALMAN  WAVE  FILTER  CALCULATIONS  AND  THEIR  REALIZATION 


Classical  Kalman  wave  filtering  is  linear  least  square 
variance  wave  filtering.  When  systems  satisfy  the  conditions 
below,  only  then  are  wave  filter  solutions  optimal. 

(1)  System  dynamics  models  and  measurement  models  are  both 
linear; 

(2)  System  dynamics  models  and  meaurement  models  match  up 
with  reality; 

(3)  Initial  state  conditions  and  a  priori  statistical 
characteristics  associated  with  noise  models  are  both 
zero  mean  value  distributions  of  already  known 
variances. 

If  the  conditions  discussed  above  are  not  satisfied,  the 
results  obtained  will  only  be  capable  of  being  approximately 
optimal.  Due  to  the  fact  that  the  majority  of  engineering 
problems  are  associated  with  dynamic  systems  and  measurement  /27 
systems  which  are,  in  and  of  themselves,  not  linear,  as  a  result, 
there  is  a  need  to  opt  for  the  use  of  generalized  Kalman  wave 
filtering. 

Assume  the  continuous  nonlinear  dynamic  equation  is 

»/<X«))  +  W,  (9) 


In  equations,  X(t)  is  the  n  dimensional  state  vector;  Wt  is  the 
m  dimensional  model  noise  vector.  It  is  zero  mean  value  Gaussian 
white  noise.  Moreover, 


E{Wt)  =  0 

E{W,  •  WJ)  =  Q(OS(t  -  t') 


(10) 


In  equations,  Q(t)  is  the  state  noise  matrix. 

Assume  again  that  the  measurement  equation  is  the  dispersion 
nonlinear  equation 

Zk+1  =  H(Xk+l)  Vk+l 
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In  equations,  Zk+1  represents  observation  vectors  at  instant  k+l. 
Vk+l  is  measurement  noise.  It  is  zero  mean  value  Gaussian  white 


noise 


E{Vk)  *  0 
E{Vk  •  V?}  =  RkSu 


(11) 


In  equations,  Rk  a  measurement  noise  matrix. 

Wave  filter  recurrence  equations  are  [7]: 

(1)  Calculation  Prediction  Values 


X1+1|1  =  X,,*  +  /(X*|i)7’  +  A(Xkik)f(Xkik)  y 


(12) 


In  this, 


A(x)  = 


3/(x) 

2r 


T  is  the  sampling 
interval . 


(2)  Calculation  Prediction  Error  Variance  Matrix 


Pk+Ut  -  <PkPk»&k  +  Qh 


(13) 


In  this,  4>k  =  I  +  A(Xk]iyr 

(3)  Calculation  Gain  Matrix 

Kk+ 1  =  Pk+utRU+i H"  R*+i) 

In  this,  Bk+x  =  ii(£»+ii») 


(14) 


if(x)  = 


3H(X) 

ax 


(4)  Calculation  Wave  Filter  Error  Variance  Array 

Fi+in+i  =  U  ~  A1+1i}i+l]Pi+ni  (15) 

(5)  Calculation  Wave  Filter  Values 

Xk+l ,J+l  =  Xk+Uk  +  X1+1[Z1+1  -  H(X1+111)]  (16) 
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Initial  recurrence  values  ^010*^010 

As  far  as  the  use  of  Kalman  wave  filter  algorithms  in  orbit 
determinations  is  concerned,  consideration  should  be  given  to  the 
several  questions  below. 

(1)  Selection  of  State  Variables 

Speaking  of  state  variable  selection  with  regard  to  Kalman 
wave  filtering  is  an  important  question.  Three  dimensional  user 
positons  and  velocities  are,  of  course,  state  variables.  Due  to 
the  fact  that  GPS  systems  make  use  of  user  clocks  in  order  to 
determine  electric  wave  arrival  times,  converting  them  into 
ranging,  and  thereby  carrying  out  positioning,  clock  differences 
and  frequency  differences  will  directly  bring  with  them  errors  in 
ranging  and  speed  determinations.  As  a  result,  user  clock 
differences  and  frequency  differences  will  also  act  as  state 
variables.  Speaking  with  regard  to  low  orbit  spacecraft,  the 
principal  perturbation  forces  influencing  their  orbits  are  the 
gravitational  forces  associated  with  the  nonspherical  shape  of 
the  earth  and  atmospheric  drag.  Moreover,  comparing  various 
other  types  of  perturbation  forces  to  these  two,  they  are  all 
high  order  and  small  in  amount.  They  can  be  regarded  completely 
as  model  noise.  As  far  as  gravitational  forces  associated  /28 
with  the  nonspherical  nature  of  the  earth  are  concerned,  it  is 
possible  to  make  use  of  mathematical  models  to  compare  and 
accurately  describe  them.  However,  atmospheric  drag,  by 
contrast,  changes  continuously  due  to  spacecraft  attitudes.  As  a 
result,  in  order  to  accurately  determine  spacecraft  positions  and 
velocities,  aerodynamic  factors  should  be  regarded  as  a  state 
variable.  In  this  way,  9  state  variables  are  included  in  wave 
filter  equations. 

X  —  \_x,y,z,x,'y,z,b,f,d'Y 

(2)  Dynamics  Models  Associated  with  Wave  Filter  Devices 

Satellite  motion  equations  are 


. 

<u  ' 

X  = 

3x 

. 

aj 

y  = 

3y 

<V 

z  = 

dz 

In  these,  U  =  (GM/r)  +  R;  GM  is  the  gravitational  constant  of 
the  earth  (3986005xl08m3/s2)  ;  r  is  the  length  from  the 
coordinate  origin  point  to  the  center  of  mass  of  the  user;  R  is 
perturbation  functions. 

In  this,  only  giving  consideration  to  gravitational  forces 
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associated  with  the  nonspherical  shape  of  the  earth  and 
atmospheric  drag,  spacecraft  motion  equations  can  then  be  written 
as 


GM  ,  ,  ,  . 

x  = - p- x  +  Au  +  Au 

y  = — —y  +  a,.  +  a2, 

GM  ...  . 
z  =■ - -p-z  +  -A  +  Au 


(18) 


% 


i 


In  this,  Alx  ,  Aly  ,  and  Alz  are  three  directional  accelerations 
given  rise  to  by  perturbations  in  gravitational  forces  associated 
with  the  earth;  A2x,  A2y  ,  and  A2z  are  three  directional 
accelerations  given  rise  to  by  atmospheric  drag  perturbations. 

(1)  Gravitational  Perturbation  Acceleration  A1  Associated 
with  the  Nonspherical  Shape  of  the  Earth 

The  potential  function  associated  with  the  earth  is  [8] 


In  this, 


U  =  Ejr)' 'fJWnW-cos-M  -  A_)  (19) 

*  - ,re,e  ZP  i  7? 


m  is  satellite  mass; 

A  is  satellite  geographical  longitude; 

ae  is  WGS-84  ellipsoid  hemi  major  axis  (6378137m) [10] ; 

Jnm  and  Anm  are  harmonic  coefficients. 

Going  through  derivations  and  appropriate  choices, 
perturbation  accelerations  associated  with  the  earth's 
gravitational  forces  are 


(20) 
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In  this,  J20  is  the  gravitational  force  field  second  order  band 
harmonic  term  10.8263x10-8  .  /29 

(2)  Atmospheric  Drag  Perturbation  Acceleration  A2 
Atmospheric  drag  perturbation  acclerations  are  expressed  as 


^PVn[_YL  =d 

2m  V'j 


v’H  - 


In  this  d  is  aerodynamic  factors; 

e  is  atmospheric  density; 

V'  is  satellite  motion  speed  relative  to  the  atmosphere. 
Going  through  resolution,  it  is  possible  to  obtain 

Au  =  —  pV’(x  +  yw,)d 
Al}  =  —  pV'Cy  —  Tw,)d 
Au  =-  pV'zd 


(22) 


The  analysis  discussed  above  is  based  on  an  execution  within 
an  inertial  coordinate  system.  However,  the  WGS-84  coordinate 
system  which  the  GPS  system  opts  for  the  use  of  is  an  earth 
fixed  coordinate  system.  As  a  result,  it  is  necessary  to  carry 
out  the  transformations 

Q 

x  —  —  pjx  4-  Au  +  Au  +  wfx  +  2w.y 
y  =  —  py  +  4-  Aty  +  w]y  —  2w.  'y  -  (23) 

z  =  —  ^z  4-  Au  +  Au 


(3)  Clock  and  Aerodynamic  Factor  Models 

On  the  basis  of  dynamics  equations — using  clock  deviation  b, 
frequency  deviation  f,  and  aerodynamic  factor  d  as  state 
variables — systems  should  satisfy 

b  —  f  4- 

+  (24) 

d  — - —d  +  tvd 


In  this,  rf  and  rd  are  first  order  Maerkefu  (phonetic)  process 
correlation  times,  wb  ,  wf  ,  and  wd  are  all  Gaussian  white  noise. 

(4)  Nine  State  Wave  Filter  Device  Dynamics  Model 
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X  =  [ x,y,z,x,y,z,b,f,dy 

x  =  /[xco]  +  w 


/30 


(5)  Precise  Determination  of  Model  Noise  Variance  Matrices 
Q 

Seen  from  the  view  of  general  references,  as  far  as  model 
noise  variance  matices  Q  are  concerned,  it  is  possible  to  select 
appropriate  values — on  the  basis  of  problems — relying  on 
experience  and  experimentation.  Among  those  commonly  used  are 
the  several  types  discussed  below  [9]: 

—  opting  for  the  use  of  rigorous  dynamics  models  in  order 
to  do  calculations  of  state  vectors; 

—  going  through  simulation  calculations  to  precisely 
determine  relatively  appropiate  Q; 

—  going  through  calculation  results  as  a  whole  to 
precisely  determine  and,  in  conjunction  with  that,  adjust 
Q,  making  differences  between  Kalman  wave  filters 
relatively  small  and  selecting  Q  for  this  time  to 
function  as  dynamic  noise  covariance; 

—  with  regard  to  maneuvering  targets,  opting  for 
the  use  of  self-adjusting  wave  filters. 

The  influences  of  dymanics  model  errors  on  orbits  are 
primarily  dependent  on  the  length  of  orbit  segmental  arcs.  The 
errors  associated  with  long  segmental  arc  orbits  primarily  come 
from  dynamics  model  uncertainties.  Speaking  in  terms  of  short 
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segmental  arcs,  the  influences  of  dynamics  model  errors  are 
relatively  small,  and  it  is  primarily  observation  noise  which 
influences  orbit  results.  As  a  result — speaking  with  regard  to 
low  orbit  spacecraft  orbits  which  are  relatively  stable — Q  can  be 
selected  as  a  constant  value. 

(3)  Precise  Determination  of  Measurement  Equations  and 
Acutalization  of  Generalized  Kalman  Wave  Filter  Devices 

GPS  goes  through  the  carrying  out  of  measurements  of  psuedo 
ranges  and  pseudo  range  speeds  with  regard  to  four  navigation 
satellites  in  order  to  calculate  user  positions  and  velocities. 
Thus,  the  observed  quantities  are  pseudo  ranges  and  pseudo  range 
velocities.  Moreover,  these  two  measured  values  are  nonlinear 
with  regard  to  state  variables  .  For  this 

reason,  they  should  be  used  to  take  a  certain  nominal  state 
expansion  and  turn  it  into  a  Taylor  series.  Omiting  higher  than 
second  order  terms,  a  linear  measurement  equation  is  constructed. 
Moreover,  nominal  states  can  be  selected  to  be  X1+1|1 
Linearization  is  carried  out  respectively  for  pseudo  ranges  and 
pseudo  range  velocities  in  regard  to  the  four  satellites,  thereby 
constructing  linearized  measurement  equations.  Having  wave 
filter  device  state  models  and  measurement  models,  it  is  then 
possible  to  realize  Kalman  wave  filtering.  In  conjunction  with 
this,  calculations  of  spacecraft  orbits  are  carried  out. 

5  SIMULATION 

This  article  opts  for  the  use  of  Monte  Carlo  methods  to 
carry  out  orbital  determination  calculation  simulations.  Results 
clearly  show  that — when  one  possesses  real  time  observation  data- 
-use  is  made  of  minimum  square  methods,  and  it  is  possible  to 
obtain  relatively  high  positioning  accuracies — or  =  15.587m,  av  = 
0.056m/s,  at  =  0.001ms.  Making  use  of  Kalman  wave  filter 
methods,  it  is  possible,  in  the  same  way,  to  achieve  relatively 
high  accuracies — ar  =  17.048m,  av  =  0.048m/s,  at  =  0.0025ms  (These 
simulation  results  do  not  consider  the  influences  of  electric 
wave  propogation  nor  do  they  consider  SA  influences) . 
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